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We propose toy models of unconventional magnetic alloys, in which the density of band states, 
p(e), and hybridization, t{e), are energy dependent; it is assumed, however, that t^(e) oc p~^(e), 
and hence an effective electron-impurity coupling r(e) = p{e)t^{e) is energy independent. In the 
renormalization group approach, the physics of the system is assumed to be governed by r(e) only 
rather than by separate forms of p(e) and t{e). However, an exact Bethe Ansatz solution of the toy 
Anderson model demonstrates a crucial role of a form of inverse band dispersion k{e). 



PACS numbers; 72.15.Qm, 75.20.Hr 



The Kondo problem in "unconventional" Fermi systems, where an effective density of states of band electrons 
vanishes either precisely at the Fermi level ( "gapless" systems) or on some interval around the Fermi level ( "gcypped" 
systems), has been attracting a significant theoretical interest. Using poor-man's scaling, Withoff and Fradkintl have 
found that the Kondo effect in gapless systems takes place only if an effective electron-impurity coupling exceeds 
some critical value. NumajiCjal renormalization group (RG) calculations, large- iV studies, and quantum Monte Carlo 
simulations of the gaplessBlJ and gappecullj systems have confirmed this prediction and revealed a number of new 
additional features of the physics of unconventional magnetic alloys. 

In a conventional metallic system with (i) a linear dispersion of band electrons near the Fermi level and (iiV-an-epergy 
independent electron-impurity coupling, ba|Si|i-Lmpurity models are exactly solved by the Bethe ansatz (BA)tj E3. It is 
well known that the Wilson's numerical RGO'cEl and BA methods lead-to identical repilts. Moreover, renormalizability 
of the Kondo and Anderson models proved respectively by AbrikosovO and HaldanenS has been reproven in the course 
of a BA solution. Therefore, it is reasonable to study the Kondo problem in unconventional Fermi systems, where the 
BA technique cannot be straightforwardly used, by making use of scaling arguments and the numerical RG approach. 

However, it has recently been founded that integrability of the U ^ oo nondegenerate and degenerate Anderson 
models is not destro5pd because of a nonlinear dispersion of particles and an epargy dependent hybridization, but it 
becomes only hiddencJ. The developed BA approach has allowed us to studycS the thermodynamics of a J7 — > oo 
Anderson impurity in a BCS superconductor, and can be used to obtain an exact solution of the Kondo problem in 
other unconventional IfefpSii systems. 

The results obtainecuJ'Ej demonstrate a discrepancy between the RG and BA approaches to the Kondo problem in 
unconventional Fermi systems. In the RG approach, the separate forms of the density of states, p(e), and an energy 
dependent hybridization, i(e), are assumed to be unimportant. It is assumed thus that the impurity properties are 
governed only by an effective particle-impurity coupling r(e) = p{e)t^{e). However, BA equations contain both an 
effective coupling and an inverse band dispersion fc(e). The latter describes the spatial behavior of wave functions 
and enters BA equations via periodic boundary conditions imposed on system's eigenfunctions. A form of inverse 
dispersion plays an extremely important role in the physics of the system. 

In this paper, we study a toy model describing an Anderson impurity embedded in a gapless Fermi system. In the toy 
model, the density of states and hybridization are energy dependent. It is assumed, however, that i^(e) = 2Tp~^{e), 
where F = const, and an effective coupling is thus energy independent. Therefore, in the framework of RG approach, 
physical properties of the toy and Anderson models are obviously identical. 

We assume a simple form of the density of states, 

= (1) 

where the energy e is taken relative to the Fermi level. The parameter P determines the size of a region with an 
unconventional behavior of the density of states. At /? = 0, we come back to the Anderson model. In the region near 
the Fermi level, |e| <C /?, the density of the states exhibits a power-law variation, p(e) ~ |e|''. 

The toy model is obviously integrable. It is clear also that a corresponding integrable toy version can be constructed 
in the described-above manner for any conventional integrable impurity model. The physical properties of the toy 
versions of the exchange models, such as the s-d (Kondo) and Cogblin-Schrieffer models, are not affected bjL-aJorm of 
dispersion, because in these models charge and spin excitations of a system are decoupled from each otheiEJ^Ej. The 
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physics of the toy versions of the Anderson models, where charge and spin excitations strongly interact, is shown to 
be governed by a form of inverse dispersion fc(e). 

We start with the Hamiltonian of a nondegenerate Anderson impurity written in terms of the Fermi operators cj. (e) 
(C(j(e)) which create (annihilate) an electron with a spin a =|, [ in an s-wave state of energy e, 



where 



H = H, + Hd + Hh (2a) 



Hc = f^^^^Ue)c4^) (2b) 

Hd = eadUa + Ud\d^dldi (2c) 
rD ^ 

H,,= &v^[4(eK+4c.(e)]. (2d) 

J-D 



are the conduction band, impurity and hybridization terms, respectively. All notation in Eqs. (2) are standard. 
The electron energies and momenta are taken relative to the Fermi values, which are set to be equal to zero. The 
integration over the energy variable e is restricted by the band half width D. In what follows, we assume that D is the 
largest parameter on an energy scale, D —^ oo. In the energy representation, the effective particle-impurity coupling 
r(e) = p{e)t'^{e) combines the density of band states, p{e) = {de{k)/dk)~^, and the energy dependent hybridization 
tie). 

Here we consider a toy version of the model assuming that i^(e) = 2Tp ^(e), where F — const. The toy model is 
obviously integrable at an arbitrary U. Indeed, let us introduce the Fourier images of the electron operators, 



c^(t) ^ J ^c<^(e) exp (ier). 



and transform thus Eqs. (2) to the auxiliary r space. In this space related to the particle energy, the model 
Hamiltonian coincides with the Hamiltonian of the standard Anderson model written in terms of the operators in the 
auxiliary x space, 

c„{x) ^ J ^c^{k)ex-p{ikx), 

related to the particle momentum. Therefore, in the t space the iV-particle wave functions of the system, 
*i'cri...(Tjv ('''ii ■ ■ ■ tTn), are given by the standard Bethe ansatz formulae derived by Wiegmanrcj. However, periodic 
boundary conditions must be imposed on a wave function '^ai...tTN (^ii ■ • ■ i ^n) on an interval of size L in the x space 
rather than in the r space. These different representations of wave functions are related by 

*o-i...o-„(a;i, • • ■ ,a;Ar) = / dri . . . tat ii(a;j- |T,)$o-^...cr„ (n, . . . , Tjv), (3a) 

where the "dressing" function u{x\t) is found to be 

1/2 

271 \ dk ) 



u{x\t)^1 ^f^V exp[^(fcx-e(fc)r)]. (3b) 



In the Anderson model, where e{k) = k, the dressing function is nothing but the Dirac delta function, u{x\t) = 5{x—t), 
and hence the x and r representations coincide. 

In what follows, we confine ourselves to the case of C/ — > oo. Then, imposing periodic boundary conditions on the 
wave function ^'cri...criv (^ii ■ ■ • tXjsi) results in the following BA equations: 

^»fc,L ^3 -ed/2F-z/2 ^ yr ujj - Aa - i/2 

ojj - ed/2r + i/2 j^j-^ujj - Xa + i/2 ^ ' 

n Xg - ujj - i/2 ^ _ TT K- Xf3-i .^l^x 
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where M is the number of particles with spin "down", kj = k{ujj), and uj = e/2T. The particle- impurity and effective 
particle-particle scattering amplitudes coincide with those in the Anderson model, because they are really determined 
only by the effective particle-impurity coupling F, as it is assumed in the RG approach. However, the BA equations 
contain also the phase factors exp {ikjL) accounting for the spatial behavior of wave functions. A form of the inverse 
dispersion k{(jj) is clear to play a very important role in both further mathematical BA constructions and in the 
physics of the system. 

As in the Anderson model, in the thermodynamic limit spin "rapidities" Xa are grouped into bound spin complexes 
of size n, 

AL"'^') = A„+z(n + l-2j)/2, ,7 = l....,n (5) 

Apart from charge excitations with real charge "rapidities" wj, the spectrum of the system contains also charge 
complexes with complex rapidities 

= A„ ± i/2, (6a) 

provided the signs of the imaginary parts of cj^^ and corresponding momenta A;^^^ = k{u)^^) are the same, 

sign(Imfc^^^) = sign(Ima;'^^). (6b) 

In the Anderson model, due to the linear dispersion law k = 2rui the necessary condition (NC) (6b) is obviously 
satisfied at an arbitrary A € (—00,00). In the toy model, a solution of NC is governed by a form of k{Lu). 

To solve Eq. (6b), one has to specify first the power r in the expression for the density of states (1). Here, we 
consider two most important cases r = 1 and r = 2. At r = 2 the inverse dispersion is found to be 

— = w — arctan — , (7a) 

where 5 = /3/2r. Solving Eq. (6b), we find a critical value of the parameter 6, 

At 6 < Scr, NC has a solution at all A. At 5 > Scr, NC has no solution on the interval Ga = {—^, where 

Thus, in a sharp contrast to the Anderson model, the spectrum of charge complexes of the system at ^ > Scr contains 
a gap of size 2 A. The gap appears a± S = Scr, and grows with increasing S. At very large S, S ^ 1, the gap 
asymptotically reaches the maximal value 2A.max — 

In a gapless Fermi system with the power r = 1, the inverse dispersion of band states is found to be 

— =w + (51n— w<0. (8) 

Solving Eq. (6b), we find no gap in the spectrum of charge complexes. Thus, the toy model with the power r = 1 
is qualitatively equivalent to the Anderson model. Therefore, in what follows we confine ourselves to the toy model 
with r = 2. 

The thermodynamics of the system is described by a set of basic cqiiations for the rcnormalizcd energies of elemen- 
tary excitation £{uj), ^(A) and K;n(A), corresponding to unpaired charge excitations with real w, charge complexes and 
spin complexes of size n, respectively: 

00 

£{uj) = 2Tu: + ai * F[-^(w)] - ^ a„ * F[-K„(a;)] (9a) 
^(A) = 4rA + ai * F[-e{\] + * F[-i{\)] (9b) 

00 

«n(A) = ^ Anm * F[-K„^{X)] + Qn * F[-e{X)]. (9c) 
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In these equations, F[f{x)\ = Tin [1 + exp (/(a;)/T)], a„(a;) = (2n/7r)(n^ + Ax^) ^, and yl„m(x) = a\n-m\{x) + 
cin+m{x) + 2^™"*-"'™^ ^ cL\n-m\+2k{x) ■ The Symbol * stands for the convolutfon of functions, e. g. 

a„ * T[C(A)] = j dX'aniX - A')FK(A')]. (10) 

For the t/ — > cx) Anderson model, the thermodynamic BA equations were derived by Schlottmann0. rlirey also can 
be obtained by setting t/ — > cx) in the thermodynamic BA equations of the general Anderson modelE3. In the toy 
model, the only difference is the appearance of the gap in the spectrum of charge complexes at S > Scr- Therefore, 
the integration contour C in the integrals with the function ^(A) consists of two intervals, C = (— oo, — A) (A, oo). 

The physical properties of the ground state of the system are governed by the gap size. At A = 0, the ground 
state, as in the Anderson model, is composed of charge complexes only. They occupy all states from A = —D/2T to 
A = Q. where Q is determined by the condition ^(Q) = 0. In conventional metallic systems Q is a large negative 
valuellj, Q ~ — (l/27r) In {D /T). The ground state of the toy model is not affected by the gap until the gap size does 
not exceed \Q\. When A exceeds Ac- = IQI, the ground state of the system is reconstructed: charge complexes with 
A € (—A, Q) decay into unbounded charge excitations and spin waves. 

Therefore, in the continuous limit the BA equations (4) describing the ground state of the system take the form 
of integral equation for the "particle", p{u}) {p{uj) = at a; > B), a{X) (cr(A) = at A > A), ?7(A), and "hole", p{uj) 
{p[ijj) = at w < B), a{\) ((t(A) = at A < —A), densities of distributions of charge excitations, charge complexes 
and spin waves, respectively, 

' -ai{u- —)^ p[Lo)+ p[uj) + ai*a(uj)+ai*ri{uj) (11a) 



2iT duj L ' 2r' 
1 ^9(A) , 1 ^^^^ - ^) = a(A) + a(A) + ai * p(A) + a2 * a{X) (lib) 



2tt dX L ' 2T' 

ai* p{X) = T]{X) + a2*r]{X), (11c) 

where q{X) — k{X + i/2) + k{X — i/2) is the momentum of the charge complexes. The "Fermi level" of unbounded 
charge excitations, B, is found from the condition 
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dojp{uj) + 2 / dXa{X). (12) 



In Eqs. (11), the densities can be separated into the host and impurity parts, e. g. p{uj) — Phi}^) + L ^pi{uj). The 
population of the impurity level, n^, and the impurity spin Sf are then given by 

Ud — I dujpi{uj) + 2 / dXai{X) (13a) 
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5f = i / diup^iiu) - I dAr,,(A), (13b) 



It is easy to see that Eqs. (11) and (13) are analogous to equations describing the ground state of the Anderson 
model in an external magnetic field Ti.. The latter, however, contain no spin waves, ri{X) — 0. Therefore, the impurity 
spin is given by 



1 f^' 



where the limit B' is determined by the field Ti. In the limit 7^ — > 0, B' ^ — oo, and the impurity spin vanishes. 

In the toy model, unbounded charge excitations and spin waves appear in the ground state of the system in the 
absence of an external magnetic field due to a decay of charge complexes. However, it is easy to show from Eqs. (11c) 
and (13b) that their contributions to the impurity spin precisely compensate each other. Thus, the Kondo effect takes 
place in the toy [/ — > oo Anderson model with the powers r = 1 and r = 2 at any effective coupling F and gap size. 

At A < IQI, unbounded charge excitations and spin waves disappear from the ground state, and the population 
of the impurity level is determined by a single function cri(A), = 2 J^^dXai{X). As in the Anderson model, the 
impurity population is governed by the renormalized impurity level energy = + 2TQ. At A > \Q\, the limit Q 
is replaced by —A, and moreover a contribution of unpaired charge excitations appears in Eq. (13a). Therefore, we 
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should expect essential changes in the behavior of the impurity population and the impurity magnetic susceptibility 
compared to the Anderson model. 

In summary, we presented an exact BA analysis of the toy version of a model describing a C/ — > oo Anderson 
impurity embedded in a gapless Fermi system. In the RG approach, the physics of the system is assumed to be 
governed by an effective electron-impurity coupling only. Therefore, the toy and Anderson models should be identical 
to each other. The BA analysis demonstrates, however, the qualitatively different behaviors of these models. 

In the Anderson model, the ground state is composed of charge complexes only. The toy model with the power 
r = 2 exhibits a critical value dcr of the parameter S describing an unconventional behavior of the density of states 
given in Eq. (1). At 5 > Scr the spectrum of charge complexes contains a gap. If the gap is quite small, A < \Q\, 
it does not effect on the ground state properties, however its appearance changes the thermodynamics of the system. 
A large gap, A > \Q\, not only changes the thermodynamic properties, but reconstructs also the ground state of the 
system, because a part of charge complexes decay into unbounded charge and spin excitations. The behavior of the 
impurity level population is drastically changed, however the Kondo effect takes place at any gap. 

Finally, it should be emphasized that the toyi-versions are interesting not only as exactly solvable examples of 
unconventional magnetic alloys. It can be shownEa that an exact BA solution of a model with the density of states 
p{e) given in Eq. (1) and an energy independent hybridization, t = const, exhibits many of qualitative features 
described above. 
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